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A 2-D impedance eduction methodology is extended to quasi-3-D sound fields in uniform or shearing mean 
flow. We introduce a nonlocal, nonreflecting boundary condition to terminate the duct and then educe the 
impedance by minimizing an objective function. The introduction of a parallel, sparse, equation solver signif- 
icantly reduces the wall clock time for educing the impedance when compared to that of the sequential band 
solver used in the 2-D methodology. The accuracy, efficiency, and robustness of the methodology is demon- 
strated using two examples. In the first example, we show that the method reproduces the known impedance 
of a ceramic tubular test liner. In the second example, we illustrate that the approach educes the impedance 
of a four-segment liner where the first, second, and fourth segments consist of a perforated face sheet bonded 
to honeycomb, and the third segment is a cut from the ceramic tubular test liner. The ability of the method 
to educe the impedances of multisegmented liners has the potential to significantly reduce the amount of time 
and cost required to determine the impedance of several uniform liners by allowing them to be placed in series 
in the test section and to educe the impedance of each segment using a single numerical experiment. Finally, 
we probe the objective function in great detail and show that it contains a single minimum. Thus, our objective 
function is ideal for use with local, inexpensive, gradient-based optimizers. 
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mode coefficients 

4MN x 1 vector of source effects 

ambient sound speed and density 

M x M matrices defining the nonlocal boundary condition 

M x 1 column vectors containing mode amplitudes and mode pressures 

M x 1 column vector containing the axial derivative of mode pressures 

d/dt + fto • V', material derivative 

transcendental function for determining eigenvalues 

source frequency 

objective function for impedance eduction 

normal derivative of 3-D source function 

3-D source functions 

height, length, and width of the duct 

x/ — 1, unit complex number 

2jt//co, freespace wavenumber 

axial propagation constant 

leading and trailing edge of test liner 

number of evenly spaced vertical and axial points 

uniform flow Mach number 

average mean flow Mach number in Grazing Incidence Tube 
3-D time dependent acoustic pressure and density 
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p,P m = 3-D and quasi-3-D time-harmonic acoustic pressure field 

p s ,p„ = source plane pressure and axial derivative of source plane pressure 

Qmn (x) = mode eigenfunction in the vertical direction 

0, ^ = normalized resistance and reactance 

[5] = 4NM x 4NM stiffness matrix 

t = time 

u, v, w = vertical, spanwise, and axial acoustic particle velocity components 
t,tf 0 = acoustic particle velocity and mean velocity vector 
wo = axial component of mean velocity 

x, y,z = Cartesian coordinates 

[Z] = dimensionless node impedance matrix 

a mn = transcendental function coefficient 

V 2 ,^ = 3-D Laplace operator and gradient operator 

Cexit = test liner and exit plane impedance 
"C; . t u = impedance of lower and upper wall downstream of the test liner 
ot = 2 ji /, angular frequency 

K„ = mode eigenvalue in the vertical direction 

{<S} = 4 MN x 1 vector of nodal acoustic pressures 

d/dn = the derivative normal to the liner surface 

• = dot product 

Superscripts: 

— 1 = matrix inverse 

* = complex conjugate 

FEM = finite element computed wall pressure 
Meas = measured wall pressure 

Subscripts: 

m, n = spanwise and vertical mode order 
0 = mean flow quantity 

J = wall location index 


I. Introduction 

E FFICIENT duct liners for broadband acoustic noise suppression are critical to the development of environmentally 
acceptable commercial aircraft. Previous experiments have shown that significant changes in duct liner impedance 
can occur as a result of grazing flow. 1 Because of the need to better understand this phenomenon, NASA Langley 
Research Center (LaRC) has invested a considerable effort into the development of finite element models for educing 
liner impedance in the presence of grazing flow. In this approach, acoustic data needed for implementation of the 
finite element methodology is acquired in the test section of the 5 1 mm x 5 1 mm Grazing Incidence Tube (GIT). This 
methodology seeks to completely describe the aeroacoustic field interior to and on the boundaries of a computational 
domain centered on the test liner (see Fig. 1). Impedances educed using this method have been extensively reported, 
and the methodology has steadily evolved from the no-flow standard of several years ago 2 to a uniform flow 3 , and 
now, to a sheared flow model. 4 To date, this methodology has been limited to 2-D aeroacoustic sound fields. 

Recently, results from the finite element method were compared to results of methods used in various U.S. aeroa- 
coustics laboratories. 5 Impedance data comparisons acquired from this multilaboratory study demonstrated the need 
to incorporate 3-D aeroacoustic effects into the impedance eduction methodologies. Subsequently, a 95-microphone 
test window was installed in the GIT 6 to obtain data for impedance eduction in 3-D aeroacoustics sound fields. Ad- 
ditionally, LaRC is building a large-scale flow duct (LSFD) that is closer in size to a full-scale engine. Because the 
cross-sectional area of the LSFD is more than 22 times that of the GIT, transverse modes (and hence 3-D acoustic 
effects) occur over the desired range of test frequencies (0.5 to 3.0 kHz). Thus, an ability to educe liner impedance in 
3-D aeroacoustic sound fields is clearly needed. 

This paper describes modifications to the LaRC 2-D finite element methodology that facilitate studies of impedance 
eduction in 3-D aeroacoustic sound fields with grazing flow. Although the method presented in this paper is applicable 
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to shearing mean flows, the results presented herein have been restricted to uniform mean flow in order to limit the 
scope of the current investigation. The paper reports initial results achieved in the GIT with these modifications using 
data acquired with a ceramic tubular test liner installed. The methodology is then applied in a geometry corresponding 
to the LSFD to demonstrate the portability of the technology to other geometries and to 3-D aeroacoustic sound fields. 
The paper also illustrates how to use the methodology to educe impedances of multisegmented liners in a flow. 

II. Statement of Problem 

F IGURE 1 shows a schematic of the 3-D duct and right-hand Cartesian coordinate system used in this study, together 
with the computational domain. The source and exit planes are located at z = 0 and z= L, respectively, as illustrated 
in the figure. The test liner is installed in the lower wall with its leading and trailing edge located at z = L i and z = Li, 
respectively. The portions of duct upstream and downstream of the leading and trailing edges are considered rigid. 
The test liner is assumed locally reacting (i.e., wave motion in the liner is restricted to the x-direction, as shown in 
Fig. 1) and has an unknown normalized time-harmonic impedance, t,, that is generally a function of wall location. 
Although the test liner is shown on the lower wall in Fig. 1, the analysis to follow will allow the liner to be placed on 
either the upper or lower wall of the duct. Throughout this work all wall impedances are normalized with respect to the 
characteristic impedance (poco) of the air flowing in the duct and all dimensions are in mks units. It is further assumed 
that the time-harmonic acoustic pressure field on the wall opposite the test liner is known. (These could be provided, 
for example, by a 95-microphone window such as that discussed in Ref. 6). Finally, there is a steady flow velocity 
vector, V o, within the duct that flows from left to right at a subsonic speed. The problem at hand is to determine the 
unknown impedance, £,, of the test liner in the presence of the viscous, compressible, steady flow field. 


III. Governing Differential Equations 


T HE governing equations for the acoustic disturbances in the duct are obtained using standard linear perturbation 
theory. The procedure and assumptions are as follows: 

1. The fluid in the duct is assumed to be an ideal gas. 

2. The steady flow is obtained from the steady solution of the Navier Stokes equations or from measurements. 

3. The acoustic equations are derived assuming that each time-dependent flow quantity may be represented as the 
sum of the known steady flow and a time dependent (acoustic) component. The sum of these two fields is 
assumed to collectively satisfy the unsteady Navier Stokes equations (turbulence is included statistically in the 
mean flow). 

4. The nonlinear acoustic equations are then obtained by subtracting off the steady flow equations from the time- 
dependent equations. 

5. We further simplify the acoustic equations by assuming the acoustic process is isentropic. Thus, viscous and 
heat conduction effects are assumed important only in establishing the steady flow field, setting up the sources 
of sound, and dissipating initial transients. 

6. The resulting nonlinear acoustic equations are then linearized about the steady flow to obtain the governing 
equations for linear acoustic disturbances in the flowing fluid. 

When written in compact vector format in Cartesian coordinates (x,y, z), the linearized acoustic equations are 7 


Dp 

Dt 


Po 


-Pof^ 

D 


r Du 


Dt 


• it ) + (it • 

o 

II 

O 

Q_ 

(1) 

= Y + 

dt 

(2) 

(t.f)Tto 

ii 

i 

<1 

(3) 

P = CoP 


(4) 


3 of 15 


American Institute of Aeronautics and Astronautics 



Equations (l)-(4) constitutes a system of five equations in five unknown acoustic variables ( p ,p, u , v, w). These equa- 
tions along with appropriate initial conditions, wall boundary condition, exit plane condition, and inflow boundary 
conditions, constitute a well-posed initial boundary value problem that can be evolved to a quasi-periodic state. Given 
the time-harmonic impedance, t, this problem can be solved to uniquely determine the time-harmonic acoustic pres- 
sure on the wall opposite the test liner. Conversely, if the time-harmonic acoustic pressure on the wall opposite the 
test liner is known (e.g., measured), the unknown impedance may be iteratively determined. The procedure would be 
to cycle through the solution in the duct for each impedance. The impedance that reproduces the known acoustic pres- 
sure opposite the test liner is the impedance of the test liner. The authors have referred to this procedure as impedance 
eduction in the duct acoustics literature. 3-6 

Henceforth, we assume that the mean flow in the duct is a fully developed parallel flow. This assumption focuses 
attention far downstream of the duct entrance, so that the initial growing shear layer has vanished, and the steady flow 
velocity vector T?q is purely axial and varies only with the two transverse coordinates; that is T?o = (0, 0,w 0 (x,y)). The 
assumption of a region of parallel steady flow at some finite distance downstream of the duct entrance is a reasonable 
approximation in many laboratory situations. 

The assumption of a parallel steady flow field allows a decoupling of the axial component of acoustic particle 
velocity (vv) from the continuity of mass flow and the two transverse momentum equations: 
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The above system (Eqs. (5)-(7)) constitutes three equations in the three unknown variables, p. u, and v, that can 
be solved to obtain the acoustic disturbances in the flow. It is convenient to combine these equations into a single 
equation involving the acoustic pressure field. The motivation for considering only an equation on acoustic pressure 
is that a pressure equation requires knowledge of only the boundary pressure disturbances to generate the appropriate 
boundary conditions. On the other hand, a solution to the above system (Eqs. (5)-(7)) will ultimately require that some 
components of the acoustic particle velocity vector be known at inflow to obtain the appropriate number of boundary 
conditions. Because the inflow boundary conditions are typically measured, and the technology for measuring acoustic 
pressure disturbances is better developed than that for measuring acoustic particle velocity, this paper focuses on the 
solution to only a single pressure equation. 

To arrive at a single equation involving the acoustic pressure field, the transverse components of acoustic particle 
velocities are eliminated from Eqs. (5)— (7). This gives a single but third order partial differential equation that governs 
the acoustic pressure disturbance in the duct 7 
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Here we are interested in only time harmonic solutions for the acoustic pressure field. These are of the form 
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Thus, for uniform flow under time harmonic conditions, Eq. (5) reduces to the convected Helmholtz equation 
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It remains to be demonstrated in a future study as to whether this derivation is sufficiently computationally robust and 
can be implemented using a practical number of acoustic pressure measurement locations. For now, it is exciting to 
note that this approach offers potential for usage in impedance eduction with shorter liners (e.g., 0.051 m long) in the 
presence of much higher mean flow velocities, which can be achieved using converging-diverging inserts in the GIT. 
Over this shorter length, it may be acceptable to assume that the mean flow is parallel sheared flow. If successful, this 
would significantly enhance our impedance eduction capability. 

IV. Acoustic Boundary Conditions 

T O obtain a unique solution to the acoustic equations (Eqs. (11) or (12)) and educe the impedance, boundary 
conditions must be prescribed. These acoustic equations require a wall boundary condition at the duct walls, a 
source boundary condition at the source plane (z, = 0), and an exit boundary condition at the exit plane of the duct 
( z = L). In addition, a normal mode analysis shows that the third order equation requires a second inflow boundary 
condition. The locally reacting, impedance boundary condition for an impenetrable surface without curvature and 
arbitrary steady flow has been presented by Myers 8 
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where the acoustic impedance is set to infinity (t, = °o) along a rigid wall portion of the duct. At the source plane the 
acoustic pressure is assumed known 

P = Ps (14) 

In addition, in a parallel shearing mean flow, we will specify the derivative of acoustic pressure normal to the inflow 
boundary as the second inflow boundary condition 



(15) 


The derivative of the acoustic pressure normal to the inflow boundary (p n ) can be obtained using the 95-microphone 
window that was discussed in a previous paper. 6 At the duct termination a variety of boundary conditions that render 
the boundary value problem on the acoustic pressure field well-posed could be implemented. When only the plane 
wave mode is cut on at the exit plane, it is convenient to specify an exit impedance condition, so that the boundary 
condition will be of the form 

d P _ ik P (16) 

dz (Mo + ^exit) 

The above exit impedance boundary condition was used in our previous 2-D analysis. 3-5 Note that the exit impedance 
boundary condition (Eq. (16)) is a boundary condition with constant coefficients that is easily implemented into our 
impedance eduction procedure. However, when higher-order modes are present at the termination, the constant coef- 
ficient exit impedance boundary condition is not appropriate. 

An alternative to the exit impedance boundary condition in the presence of higher-order modes is the nonreflecting 
boundary condition. Physically, implementing a nonreflecting boundary condition at the exit plane is equivalent to 
attaching an infinite length of duct downstream of the exit boundary so that downstream propagating waves are not 
reflected. Such a condition is generally accepted as being physically correct when the portion of duct downstream of 
the trailing edge of the test liner is extremely long. In this work we apply a generalization of a zero-flow, nonreflecting 
boundary condition presented in an earlier work. 9 A detailed derivation of the nonreflecting exit plane boundary 
condition that includes mean flow effects is presented in a subsequent section. 


V. Quasi-3-D Analysis 

I N this paper it is convenient to develop the impedance eduction theory in Quasi-3-D. The quasi-3-D theory greatly 
reduces the computational time and memory required to educe the impedance by reducing the 3-D problem to a 
series of 2-D calculations. However, the quasi-3-D model still adequately represents many of the essential features of 
a fully 3-D theory. The quasi-3-D theory is applicable when the solution in the spanwise direction is separable (e.g., 
when one or both sidewalls of a duct are lined with uniform acoustic treatment and the mean flow is independent of 
the spanwise coordinate). When the sidewalls of the duct are rigid as in Fig. 1, the acoustic pressure field and sound 
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source may be expanded into the following Fourier series of hard wall duct modes 
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where m takes on integer values, and the functions G m (x) and g m {x) are known from the source condition. 

Thus, the partial differential equation, wall, and source plane boundary conditions that govern impedance eduction 
in uniform or parallel shearing flow are now summarized. In a uniform flow the differential equation is 
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and inflow boundary condition is 

Pm — G m {x) 

On the other hand, when there is shear flow, the differential equation is 
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and the inflow boundary conditions become 
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Equations (20)-(22) for uniform flow or equations (23)-(26) for parallel shearing flow, along with the appropriate 
nonreflecting, exit plane boundary condition, constitute a well-posed boundary value problem for the unknown pres- 
sure function P, n {z,x). Exact analytical solutions for P m (z,x) are available for infinitely long ducts. This infinite duct 
analytical solution is derived in the following subsection and its derivation is essential for formulating the proper 
mathematical form of the nonreflecting boundary condition at the duct exit. 


VI. Nonreflecting Boundary Condition 

A Nonreflecting boundary condition at the exit boundary is now developed. This boundary condition can be used 
provided that the duct termination is nonreflecting. The nonreflecting boundary condition developed here is a gen- 
eralization of the zero flow boundary condition developed in an earlier work. 9 The development begins by expanding 
P,„ and its axial derivative at the exit plane as a series of right-moving duct modes 
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We now divide the vertical coordinate, x, into M evenly spaced nodes, located at x = x\ , x = X 2 , . . .x = xm, across the 
exit of the duct. Evaluating the series in Eqs. (27) and (28) at each boundary node point gives 
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The solution to Eqs. (29) and (30) is 

-ik{P} = [Z]{U} 
[Z] = [C][G]- 1 [C]- 1 


(31) 


(32) 


(33) 


(34) 

(35) 


Equation (34) is a nonlocal, nonreflecting boundary condition that allows all waves to exit the exit plane without 
reflecting. This boundary condition has the disadvantage that the matrix [Z] is a dense matrix that must be numerically 
computed. On the other hand, this most general form of the boundary condition contains each and every propagating 
mode, is not dependent on the sound source, and is valid in flows with or without mean shear. Note that when a single 
mode, (mn), exists at the termination, the boundary condition reduces to the well-known, local boundary condition for 
a single outgoing wave in a duct 

{U} = -iK mn {P} (36) 

The dimensionless node impedance matrix, [Z], can be fully determined once the mode shape Q m „(x) in the vertical 
direction and axial wave numbers K mn are determined. We now derive expressions for these quantities assuming that 
the upper and lower walls of the duct downstream of the trailing edge of the test liner are lined with uniform acoustic 
material. The analysis will assume uniform flow for the sake of simplicity, but it is easily extendible to the case of 
a parallel shearing mean flow. In the case of uniform mean flow, the eigenmode in the vertical direction (Q mn (x)) is 
obtainable using separation of variables 
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Note that for a given value of K„, the axial propagation constant K mn , satisfying Eq. (40), has two roots. Here, we will 
choose the root that corresponds to right moving waves in the duct. 

The method of analysis described above requires a calculation of the mode eigenvalues, that are the complex 
zeroes of the transcendental function, E(k,H ,M qXiXu^Ki)- The method used to obtain these roots parallels the 
analytic continuation method discussed in Ref. 10 and is not discussed in further detail here. It should be noted that 
in most flow impedance tubes, the section of duct downstream of the trailing edge of the test liner contains rigid walls 
(i.e., 'C„ = X>1 = 0 °) so that the zeroes for the transcendental function, E(k,H, Mo, oo, oo, are known 

tm 

K= — ,n= 0,1,2,... oo (41) 

ri 

and the axial propagation constants can be determined exactly 



VII. Impedance Eduction Technique 

T HE numerical method chosen to solve Eqs. (20)-(22) for uniform flow and Eqs. (23)-(26) for parallel shearing 
flows, along with the nonreflecting boundary condition, closely parallels that used in an earlier paper 3 for the 
zero-order spanwise mode (i.e., m = 0). It is assumed that there are N and M evenly spaced nodes in the axial and 
vertical directions of the duct, respectively. A finite element method with cubic Hermite polynomial basis functions 
and a rectangular element 11 is applied to obtain the solution for P m (z,x). Integrating the second derivative terms 
in the partial differential equation by parts allows the wall impedance and nonreflecting boundary condition to be 
incorporated into the methodology using a weak formulation. The finite element method results in a discrete system 
of equations of the form 

[5]{T>}={B} (43) 

where [5] is the sparse, complex, asymmetric stiffness matrix with order 4 MN. We use an asymmetrical, parallel, direct 
sparse solver to obtain {<!>}. The solver employs a compressed column storage scheme to reduce storage overhead. 
Only the nonzero coefficients within [5] are stored, along with two pointer arrays that store the column numbers and 
starting indices of these nonzero coefficients. To achieve efficient solutions, the sparse solver uses two accelerators: 
equation reordering 12,13 to reduce fill during the factorization of [5] and parallelization. (The equation solver runs 
on multiple processors simultaneously.) Stewart’s adaptation of the Davidon-Fletcher-Powell (SDFP) optimization 
algorithm 14 is used to determine the unknown impedance, ’C, by minimizing the objective function 

F&) = f [F^ M (zj,H)-P™™(zj,H)] [F^ m (zj,H)-P™™(zj,H)]* (44) 

. 7=1 

VIII. Results and Discussion 

A n in-house computer code (FEMQ3D) that implements the impedance eduction methodology discussed in this 
paper has been developed. The code uses F90 Fortran with double precision (i.e., 64 bit) arithmetic and was de- 
signed to run in a multiprocessor environment using a shared memory implementation. This approach has been chosen 
because the primary hardware to be utilized was a Silicon Graphics Origin 3000 system. In the current implementation 
the resistance, 0, and reactance, of the test liner are the design variables ('C = 0 + {%). The optimizer (i.e., SDFP) 
runs only in sequential mode using central finite differences to compute the gradient of the objective function. The 
finite difference step size is 1 x 10 -8 , and a stopping criteria of 1 x 10 -8 (F(£) < 1 x 10 -8 ) is used to terminate the 
search. Finally, we have accelerated the impedance eduction methodology by running the equation solver in parallel 
on eight central processing units (CPUs). 

Two examples are provided to demonstrate the accuracy, efficiency, and robustness of the impedance eduction 
methodology presented in this paper. In the first example we demonstrate the ability of the methodology to converge 
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to the known impedance of a ceramic tubular (CT57) uniform test liner. This example contains two design variables, 
the resistance and reactance of the uniform liner. In the second example, we show that the method can educe the 
known impedances of a four-segment liner, in which three segments contain conventional liner configurations (three 
different perforated facesheet-over-honeycomb core combinations) and one segment is a CT57 liner. This second 
example contains eight design variables, the resistance and reactance of each of the four liner segments. 

We desire to implement this quasi-3-D methodology for the LSFD currently being brought on-line at NASA LaRC. 
Since this duct is not yet fully operational, measured data are not available. Therefore, for the current study, we use 
synthesized data based on impedances educed from data acquired in the GIT. For the first example, the impedance of 
the CT57 liner is used as input to the exact modal solution for outgoing waves in an infinite duct (see the discussion 
at the end of section VI). To simulate the infinite duct, we set Li = 0.0 and Li = L (see Fig. 1). It should be noted that 
for each spanwise mode order (m) and frequency (/), the mode with the lowest decay rate is used as the sound source 
for this example. 

Unfortunately, exact solutions are not available to synthesize acoustic pressure data for a four-segment liner in the 
LSFD geometry. Instead, a finite element code was used to synthesize data for the four-segment liner by incorporating 
the known impedances of each of the four segments. For this study, the first, second and fourth segments were 
conventional liner segments, and the third segment was a CT57 liner segment. Each segment is four inches long 
and eight-inch sections of rigid wall are located upstream and downstream of the leading and trailing edge of the 
four-segment liner. We used a plane wave as the sound source in the hard wall duct for this example. 

Both examples (i.e., the uniform CT57 liner and four-segment liner) use a nonlocal, nonreflecting boundary con- 
dition (Eq. (34)) to terminate the duct. Also, the SDFP optimization is initialized using normalized resistance and 
reactance values of 1.0 and -1.0, respectively. Further, a uniform spatial grid with 65 points ( N = 65) in the axial 
and 17 points ( M = 17) in the vertical directions, respectively, was used to educe the impedances. Thus, the node 
impedance matrix, [Z], used in the computation was a 17 x 17 dense matrix. It should be noted that since there are 
slightly more than seven wavelengths in the axial direction (for a hard-wall duct), the chosen spatial grid is approxi- 
mately three times finer than is needed in that direction. 

The impedances for each example are educed under conditions for which the liner is located on the upper wall and 
higher-order spanwise modes are cut on. These conditions are implemented for a range of frequencies and flow Mach 
numbers. The frequency ranges from 0.5 to 3.0 kHz (in increments of 0.1 kHz) and the centerline Mach numbers are 
0.0, 0.2, and 0.4. To limit the scope of our investigation, we present results only for uniform mean flow. It is our 
intention to present results with shear at a later time. 

Before proceeding to the impedance eduction studies, it should be noted that a typical impedance eduction run (i.e., 
for a single Mach number, frequency, and spanwise mode number) using FEMQ3D on a single CPU requires slightly 
less than 40 seconds of wall clock time and 200 megabytes of RAM on our Silicon Graphics ORIGIN 3000 computer 
platform. Additionally, the new sparse solver is nearly an order of magnitude faster than the now antiquated band 
solver used in our earlier papers. 2-4 Additionally, results presented in the current paper were achieved by running 
the sparse equation solver in parallel on 8 CPUs. This reduced the wall clock time from 40 seconds (on a single 
CPU) to slightly less than six seconds (using 8 CPUs) for each optimization run. These large reductions in wall clock 
turnaround times should allow the current method to be applied to a fully 3-D problem. 

A. CT57 Liner 

T HE resistance and reactance educed in the CT57 liner (from synthesized data) in the LSFD at zero Mach number 
(Mo = 0.0) are plotted in Fig. 2. We use a dual-axis system with the resistance and reactance plotted on the yl and 
y2 axis, respectively. Results obtained for each spanwise mode, m , in the LSFD are plotted against those educed in the 
GIT. To the level of resolution presented in Fig 2, the impedance educed in the LSFD is the same for each spanwise 
mode and reproduces exactly that educed (from measured data) in the GIT. For this reason, we have plotted the GIT 
and the LSFD results for each spanwise mode order as a single curve in Fig. 2. This excellent agreement is evidenced 
by the coalescence of four curves (one for the GIT and three for the individual spanwise modes in the LSFD) for 
each component of impedance, and results from the fact that an extremely fine finite element spatial grid was used. 
It is also clear that the small values of the finite difference step sizes and stopping criterion used in SDFP are also 
important parameters contributing to such excellent agreement. The number of iterations that SDFP requires to educe 
the impedance is an increasing function of the spanwise mode order. As an example, the second order spanwise mode 
(m = 2) typically requires four times the wall clock time as that of the zero order spanwise mode (m = 0). 

Results are given in Figs. 3 and 4 for mean flow average Mach numbers (averaged across the duct cross-section) 
of 0.172 and 0.335, respectively. Recall that the impedance of the CT57 liner was obtained in the GIT where bound- 
ary layers were present. Note, this uniform flow calculation ignores boundary layer effects. These uniform Mach 
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numbers were used with the LSFD synthesized data to educe the impedance. The key feature of this exercise was 
to demonstrate self-consistency of FEMQ3D when higher-order modes are included. It is also an indication that the 
nonlocal, nonreflecting boundary condition (Eq. (34)) is functioning properly in the presence of grazing flow. Again, 
it was observed that the SDFP required many more iterations (and thus, significantly more wall clock time) to educe 
the impedance when higher-order spanwise modes were propagating in the duct. 

A major drawback introduced by the use of local, gradient-based optimizers (such as SDFP) is that they often 
converge to a local minimum. In this situation, different design points (i.e., impedances) may be obtained if the 
optimizer is initialized from different points in the design space. The reduced wall clock turnaround time of FEMQ3D 
has allowed the authors to perform exhaustive studies of the effects of the starting location. The SDFP was initialized 
at 100 distinct starting locations and was found to converge to the same value of impedance for 95 of these starting 
locations. The SDFP did not converge for the remaining five starting locations because of a tendency of the gradient 
search technique to stall on flat surfaces. This clearly indicates that the objective function (Eq. (44)) contains a single 
minimum point for the uniform liner. 

B. Four-Segment Liner 

T HE second example presented is for the four-segment liner. The current quasi-3-D methodology allows the 
impedance of each of the four segments to be educed in the presence of flow. This is a significant enhance- 
ment over the previous impedance eduction methodologies implemented by the authors. Educed impedances were 
obtained by initializing each of the four impedance segments with 'Q = 1.0— 1.0/ in the SDFP optimization. Because 
the impedances of the three conventional liner segments were only known for average mean flow Mach numbers of 
0.172 and 0.335, and for frequencies of 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 kHz, computed results are restricted to these 
ranges of parameters. 

Figures 5a, 5b, 5c, and 5d show the educed resistance and reactance values for the first, second, third, and fourth 
segments, respectively. The average mean flow Mach number is 0.172, and results are plotted only for the zero order 
spanwise mode (m = 0) and first-order (m = 1) spanwise modes. With the exception of the anti-resonance frequency 
(2.0 kHz) results for the third segment (CT57), the educed resistance and reactance of each segment is nearly the same 
(as expected) for each spanwise mode (m = 0 and m = 1). Also, these values are well matched to those educed in 
the GIT. The discrepancy in the impedance of the third segment (see Fig. 5c) at the antiresonance frequency has been 
observed in the 2-D analysis in an earlier paper 3 and is believed to be due to the sensitivity of the sound field to small 
changes in impedance near the antiresonant frequency. Segmented liner impedance eductions for an average mean 
flow Mach number of 0.335 are shown in Figs. 6a-6d. Trends at this higher Mach number are consistent with those in 
Figs. 5a-5d. 

One complication in the four-segment liner impedance eduction is the increase in the number of SDFP stall points 
when compared to that of the uniform CT57 liner. This result is not that surprising, given that the function space has 
increased from a two-dimensional to an eight-dimensional design space. We have also experimented with different 
starting locations in SDFP for the four-segment liner. Each segment of the four-segment liner was initialized from 
the following six starting locations, 1 + /, 1 — /, 5 + 5/, 5 — 5 /, 10 + 10/, and 10 — 10/. We observed that the SDFP 
converged to the same value of impedance as plotted in Figs. 5 and 6, except for the instances where the SDFP stalled 
on a flat surface and had to be restarted elsewhere in the design space. 

IX. Conclusions 

T HE LaRC impedance eduction methodology has been extended to account for quasi-3-D sound fields in a uniform 
or shearing flow. A key element of the shear flow methodology is the introduction of a third orderpressure equation 
that does not require a measurement of the inflow acoustic particle velocity for its implementation. A second important 
element of the current theory is the formulation of a nonlocal, nonreflecting boundary condition to terminate the duct 
when higher-order spanwise modes are cut on. The introduction of a parallel, sparse, equation solver has significantly 
reduced the wall clock time for educing the impedance when compared to that of the band solver used in earlier 
implementations. Results presented here have shown that the uniform flow methodology consistently reproduces 
the same impedance as that achieved for the zero order spanwise mode when several higher order spanwise modes 
(m = 1,2) are cut on, thus demonstrating inter-code consistency. This work also illustrates that the uniform flow 
methodology is capable of educing the impedance of a multisegmented liner. This ability should significantly reduce 
the amount of time and cost required to determine the impedance of several uniform liners by allowing them to be 
placed in series in the test section, and the impedance of each segment educed using a single experiment. Finally, we 
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have demonstrated the robustness of using the known acoustic pressure on a wall opposite the test liner to construct 
the objective function for the impedance eduction. This objective function contains only a global minimum and is, 
therefore, ideal for use with local, inexpensive, gradient-based optimizers. Also, the sheared flow impedance eduction 
methodology needs to be tested and validated. 
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Figure 1. Schematic of 3-D duct, computational domain, and acoustic test liner. 
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Figure 2. Educed impedance for CT57 liner in LSFD with synthesized data compared to measured data in GIT (Mo = 0.0). 
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Figure 3. Educed impedance for CT57 liner in LSFD with synthesized data compared to measured data in GIT (Mo = 0.172). 
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Figure 4. Educed impedance for CT57 liner in LSFD with synthesized data compared to measured data in the GIT (Mo = 0.335). 
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Figure 5. Educed impedance of four-segment liner in LSFD with synthesized data compared to measured data in GIT at Mach 0.172. 
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Figure 6. Educed impedance of four-segment liner in LSFD with synthesized data compared to measured data in GIT at Mach 0.335. 
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